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ABSTRACT 

We describe how thermally produced anti-electron neutrinos (P e ) from the 
homologously collapsing core of a supermassive star (M ^ 5 x 1O 4 M ) can 
lead to significant deuterium enrichment in the ejected envelope of such a 
star. Deuterium-enriched material at high redshift might then serve as a clue 
to the existence of pregalactic supermassive stars. Conceivably, the ejected 
deuterium-enriched material could intercept the line of sight to a distant QSO 
and mimic a Lyman limit absorber. In such a case, the deuterium abundance 
inferred from absorption lines might not reflect the true primordial abundance 
of deuterium. We discuss relevant theoretical uncertainties in supermassive star 
physics as well as potential observational signatures in Lyman a absorber clouds 
for processing by stars of this kind. 

Subject headings: elementary particles - nuclear reactions, nucleosynthesis, 
abundances - cosmology: observations and theory 
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1. Introduction 

In this paper we assess the prospects for deuterium synthesis resulting from z/ e -induced 
neutron production in supermassive stars (M ^ 5 x 10 4 M Q - these objects become 
dynamically unstable due to the general relativistic instability). We will show that this 
process cannot account for the observed interstellar or solar system deuterium abundance. 
However, it might result in local enrichment of deuterium above its primordial value in the 
ejecta from a collapsed supermassive star. This might give a telltale sign of the existence 
and evolution history of pre-galactic supermassive stars - as yet there is no direct evidence 
for these objects ever having been extant in the universe, though it has been argued that 
their formation may be an inevitable consequence of the evolution of dense star clusters 
and perhaps ultra-low metallicity environments in the early universe (Hoyle & Fowler 1963; 
Begelman & Rees 1978; Fuller, Woosley, & Weaver 1986, hereafter FWW; Bond, Arnett, & 
Carr 1984; McLaughlin & Fuller 1996). 

If such deuterium-enriched ejecta at high redshift intercept a line of sight to a more 
distant QSO, then it could mimic a Lyman limit absorber. Measurement of isotope-shifted 
hydrogen absorption lines in such a system may then lead to an erroneous estimate of the 
primordial deuterium abundance (D/H) and, hence, an incorrect inference of the fraction 
of the closure density contributed by baryons, Qb- Future indisputable observations of 
significantly discordant D/H values among different Lyman limit systems, or very high 
values of this quantity, may be a signal for the existence of pregalactic supermassive stars. 
Confirmation of such a conjecture would depend on other evidence of supermassive star 
processing (FWW): nucleosynthesis products of hot hydrogen burning (the rp-process, 
Wallace & Woosley 1981); greatly enhanced helium abundances; or effects of large black 
holes. 

Deuterium is very fragile in stars. Thus to search for the primordial value of D/H one 
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has to look for gas clouds unprocessed by stars. Low metallicity Lyman a absorbers at 
high redshifts are thought to be such primordial clouds with little stellar processing (cf. 
Jedamzik & Fuller 1997). So far, measurements derived from three Lyman limit absorption 
systems are available. One indicates a high value with D/H= 1.9(±0.5) x 1CT 4 (Songaila 
et al. 1994; Ruger & Hogan 1996; see however Tytler, Buries & Kirkman 1996), and two 
yields low values with D/H= 2.3(±0.3 ± 0.3) x 1(T 5 and 2.5(±0.5 ± 0.4) x 10~ 5 , respectively 
(Buries & Tytler 1996; Tytler, Fan & Buries 1996). The implications of these measurements 
for primordial nucleosynthesis and the value of Qb have been discussed extensively (Cardall 
& Fuller 1996; Hata et al. 1997; Jedamzik & Fuller 1997; Copi, Olive, & Schramm 1996). 
One may be left with an uneasy feeling regarding how representative the three systems are 
of Lyman limit systems in general. By necessity only systems with very narrow lines can 
show the sought after isotope shift, and only systems with simple line profiles can minimize 
the chance of blending. Clearly these systems comprise a small fraction of the Lyman limit 
systems along the lines of sight to bright QSO's. 

The great bulk of deuterium in the universe is generated by freeze-out from nuclear 
statistical equilibrium at high entropy in big bang nucleosynthesis (Wagoner, Fowler, & 
Hoyle 1967). It has been argued effectively that the bulk of the deuterium cannot be 
synthesized through cosmic ray spallation or photo-disintegration of nuclei (e.g., Epstein, 
Lattimer & Schramm 1976; Sigl et al. 1995). 

There are two characteristics of the deuteron which make its synthesis problematic: 
(1) it is the most fragile nucleus, with a binding energy of only 2.2 MeV; (2) it requires 
free neutrons. It is not easy to obtain free neutrons in astrophysical environments, and 
to complicate matters, neutrons are unstable. They can only be "spalled" from bound 
nuclei, or produced directly via the weak interaction. The former production channel is 
difficult because neutrons are bound in nuclei by ~ 8 MeV; while the latter channel is 
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problematic because weak interactions are inherently slow. Having somehow produced a 
free neutron, significant synthesis of deuterons then requires that the n(p,7)D rate be fast 
compared to the free neutron decay rate and the local material/thermodynamic expansion 
timescale, while the primary deuteron destruction process, D(p,7) 3 He, is slow compared to 
the local material/thermodynamic expansion timescale. Supermassive stars come close to 
satisfying these requirements as they are low baryon density, relatively cool, high entropy 
configurations, in which the n(p,7)D rate is relatively fast and the D(p,7) 3 He rate is 
relatively slow. The free neutrons, in this case, are produced by intense neutrino fluxes 
released during the collapse of supermassive stars via the weak reaction v e + p — > n + e + . 

The crux issue in this scenario for deuterium production hinges on whether any 
material in these stars can be ejected from regions where all the requirements to yield 
significant D/H enrichment are satisfied. Indeed, Woosley (1977) suggested that something 
like this process could account for all of the deuterium in the universe. Though we will 
argue that this is not possible given our current understanding of how supermassive stars 
evolve (FWW), it may still lead to significant local deuterium abundance enrichments, to 
as high as D/H ^ 1CT 4 . We also show that the possibility that supermassive star ejecta 
form a "Lyman limit system" is consistent with the characteristics of the currently observed 
Lyman a absorption systems from which D/H values have been measured. We then discuss 
how to eliminate this possibility by further observations. 

2. Deuterium Enrichment in the Envelopes of Supermassive Stars 

In this section we will describe how D/H can be enhanced in some regions of a 
supermassive star via this three step scenario: 



• The star collapses and, near the point where the homologously collapsing core 
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becomes a black hole, a fraction of the gravitational binding energy is radiated as 
thermally-produced vv pairs. 

• The process v e + p — > n + e + creates free neutrons further out in regions of the star 
where material ejection is possible. 

• These neutrons undergo n(p,7)D on a timescale short compared to the material 
expansion timescale and the neutron decay lifetime; yet the lifetime of the deuteron 
against destruction via D(p,7) 3 He is long compared to the disassembly time (or 
inverse expansion rate) of this region of the star. 

Leaving aside the thorny and unresolved issue of how a supermassive star could 
be formed, it is instructive to review the simple structure and evolution of nonrotating, 
primordial-composition supermassive stars (FWW). These stars are radiation dominated, 
with essentially all the pressure support coming from photons, while nearly all the 
mass-energy in the star is contributed by baryon rest mass. These stars will have an 
adiabatic index Y\ 4/3, will be nearly constant entropy, fully convective configurations, 
and as a result will be very accurately described as an index n = 3 Newtonian polytrope. 
Such a configuration will have close to zero total (internal plus Newtonian gravitational) 
binding energy and will be "trembling on the verge of instability" as outlined by Fowler 
(1964) and Iben (1963). Additionally, this configuration will be radiating at the photon 
Eddington limit. As a result, the star will resemble an over-sized Wolf-Rayet star and will 
be accompanied by continual and significant mass loss during its lifetime. This lifetime will 
be quite short, however, as the star will quasi-statically contract as it radiates away entropy. 
This will continue until the central density reaches a critical value, beyond which general 
relativistic effects ensure that the star becomes dynamically unstable. This instability point 
is roughly coincident with the onset of hydrogen burning for stars with masses ~ 10 5 M Q . 
The subsequent fate of the star depends on a number of factors, principally composition, 
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rotation and magnetic pressure. Explosion or collapse to a black hole depends on how 
rapidly thermal energy can be added to the system through the competition between hot 
hydrogen burning on the one hand, and thermal neutrino losses and the build-up of infall 
kinetic energy on the other. In general, collapse to a black hole is probably inevitable for a 
zero metallicity, nonrotating initial configuration (FWW). 

However, there is a generic feature of this collapse which serves to isolate an outer 
region of the star from a "homologous core" which will plunge through the event horizon 
as a unit. Initially the entire star will collapse homologously, but as pressure is relatively 
reduced in the center by thermal neutrino emission, only an inner part of the star can 
continue to collapse homologously (Goldreich & Weber 1980; FWW). The boundary of this 
homologous core (as well as the mass enclosed) depends on the entropy per baryon, which 
in turn depends on the history of nuclear burning-driven entropy increase and neutrino 
loss-driven entropy decrease (FWW; Fuller & Shi 1997). Material outside of this core 
could be subject to ejection resulting from neutrino heating (Fuller & Shi 1997), and/or 
hot rp-process nuclear burning-driven convection combined with rotation and MHD effects 
(c/., Ozernoi & Usov 1971; Fowler 1966). Our work here is meant to suggest that detailed 
multi-dimensional numerical modeling of these scenarios merits serious attention. 

Prodigious thermal neutrino losses will accompany the collapse of the 
homologous core (FWW). The peak neutrino-pair production rate will come 
near the point where the homologous core of mass M HC approaches its 
Schwarzschild radius, r s ps 3 x 10 cm (m hc /m ), where the average density is 



p s a* 1.8 x 10 16 gcm 3 (M HC /M ) . Crudely, the Newtonian gravitational binding energy 




at this point is Eq 



0.5M HC c 2 ~ 10 54 ergs (m hc /M ), while the free fall timescale is 
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Not all of the gravitational binding energy will be carried away by neutrinos. A 
fraction of the binding energy could reside in infall kinetic energy and thermal energy, or 
be dissipated in magnetic modes, or be "stored" in rotational kinetic energy. In addition, 
neutrinos radiated from the core will be gravitationally redshifted. We designate the 
fraction of the binding energy radiated by neutrinos by f E . This quantity /e will be a 
significant fraction of unity if neutrinos are not trapped in the homologous core as a result 
of their mean free path becoming short compared to the core radius. (Fuller & Shi 1997 
find that such trapping is unlikely except very near the black hole formation point for stars 
with initial masses above M ^ 10 5 M .) 

During the collapse the temperature will be T 9 = T/10 9 K « 0.67(5.5/ g) 1/s S^p 1 ^ , 
where S100 is entropy per baryon in units of 10 2 Boltzmann's constant, P3 is the density 
in 10 3 gcm~ 3 and g is the statistical weight in relativistic particles (g ~ 2 and g ~ 11/2 
for T 9 ^ 1 and T 9 ^ 1, respectively). Si 00 is typically between 1 to 10. Near the endpoint 
of collapse, the core temperature will be of order T 9 HC ~ 11.7(M^ C ) ^ where Mf c is 
the homologous core mass in units of 10 5 M at the blackhole formation point. Allowing 
for entropy loss by thermal neutrino emission and including a fair amount of angular 
momentum/centrifugal support for M 5 w 10 suggests that Mf c could be about an order of 
magnitude smaller than the initial stellar mass (Fuller & Shi 1997). This implies that typical 
temperatures will be T 9 ~ 10 to 20 at the collapse endpoint for initial stellar configurations 
with 1 ^ M 5 ^ 10. For T 9 ~ 10 to 20 the average energy of neutrinos produced as pairs in 
e ± annihilation will be (cf. Schinder et al. 1987; Itoh et al. 1989) (E„) ~ 6 to 12 MeV. 
The neutrino energy will suffer an undetermined gravitational redshift of order 1, but for 
simplicity we absorb the redshift into the factor / E . 
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The probability that the process v e + p — > n + e + converts protons at radius r into 
neutrons is: 

E™ ax fW 1 \/ n,, \. .. 

1 



The third term in this equation is the fraction of neutrinos which are P e , and is typically 
about 1/3. The cross section for v e absorption on protons is a {E v ) pa 10~ 43 cm 2 (i?,,/MeV) 2 
when E u > 1 MeV. Therefore, 

With this expression we can estimate that P pa 1.5 x 1CT 5 when / E pa 0.1 and r = 3 x 10 12 cm 
for a homologous core mass of Mf c = 1. This is an interesting figure of merit, as it 
is roughly the lowest reasonable value for the primordial deuterium abundance - any 
production above this value would correspond to D/H enhancement. A P of 1CT 4 is not 
out of reach if f E can be as high as 60%. Note that a baryon at r = 3 x 10 12 cm for the 
indicated value of Mf c would have a gravitational binding energy of order 0.5% of its rest 
mass. The key question is whether there are any circumstances under which material with 
this binding energy could acquire enough thermal energy that ejection would be possible 
(Fuller & Shi 1997). 

Here P will also be the amount of deuterium enrichment relative to hydrogen in the 
limit where the reaction n(p,7)D is fast compared to the material ejection rate and the free 
neutron decay rate and where the deuteron destruction process D(p,7) 3 He is slow compared 
to the material ejection rate. In fact such conditions can be attained for some outer regions 
of the supermassive star as we illustrate in Figure 1. In this figure the material ejection 
timescale is taken to be the free-fall timescale, as they are typically of the same order. The 
stippled region corresponds to thermodynamic and expansion conditions where deuterium 
can be produced and can survive. 
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The fate of the initial primordial deuterium in the envelope is unclear. Supermassive 
stars are convective during their pre-collapse contraction phase. If their convection proceeds 
faster than the contraction, and their central temperatures during most of the contraction 
phase are high enough to burn deuterium, the primordial deuterium may be destroyed in 
the entire star before collapse. This could occur in less massive supermassive stars whose 
central temperatures are higher. 

3. Envelopes of Supermassive Stars as Lyman a Absorbers 

Once the enriched envelope of a supermassive star is ejected, it expands freely in space 
and eventually cools down to ~ 10 4 K, in equilibrium with ambient UV radiation. Its 
expansion velocity is at least ~ 10 km/sec, the thermal speed at 10 4 K. It can easily expand 
to a size of 1 kpc, unless it encounters another cosmic structure or is slowed down simply 
by loading of the ambient medium. The ambient baryonic mass swept up by the envelope is 
on average (4ir/3)R 3 ilbP ~ 10 3 fi& h 2 (1 + z) 3 (i?/lkpc) 3 M , where p is the average matter 
density at a redshift z, and R is the radius of the expanding envelope. At z ~ 2 to 3, 
this mass will be larger than the mass of the ejected envelope if R ^> lkpc. Therefore, 
the expansion of the envelope up to R ~ 1 kpc can be essentially free, and its chemical 
composition can be essentially unperturbed. 

The number of such envelope systems in one quasar line of sight is 

(Q sms p/M)nR 2 l(t/t ) ~ 300a ms /i 2 (l + ^) 3 (i?/lkpc) 2 M 5 - 1 , (3) 

where f2 sms is the density parameter of baryons in supermassive stars, and M their masses. 
/ ~ 10 3 Mpc is the proper distance interval under investigation in a quasar line of sight, 
t ~ 10 8 years is the duration of the free expansion, and to ~ few x 10 9 years is the age of 
the universe then. Therefore, even if we assume that 10~ 5 of the mass of our universe (or 
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0.01% to 0.1% of baryons) is in these supermassive stars at z ~ 2 to 3, we can expect to see 
at least one such envelope absorber in tens of quasars. In comparison, there are reported 
detections from 3 Lyman limit systems from surveys of ~ 100 quasars, a statistic that can 
easily be accommodated by these absorbers being ejected envelopes of supermassive stars. 

The column densities of these ejected ejected envelopes are approximately 

xU cm /TxR 2 m p ~ lO 19 x(i?/lkpc)- 2 (M cnv /lO 5 M ) cm" 2 , (4) 

where x is the ionization fraction of the envelope, M env is its mass, and m p is the proton 
mass. With a temperature of ~ 10 4 K and the envelope in photoionization equilibrium 
with the ambient photoionizing UV background, x ~ 10~ 2 . Therefore the resultant column 
density is that of a Lyman limit system, ~ 10 17 cm -2 , with M env being a significant portion 
of the parent mass ^ 1O 5 M . 

Curiously enough, all three Lyman limit systems with reported detections of D/H have 
double components. This is fully consistent with the geometry of ejected envelopes. The 
velocity separations between the double components range from 11 to 20 km/sec, which can 
easily be produced if the envelopes expand at ~ 10 km/sec. The double components that 
yield the lower D/H are quite asymmetric, however. But one cannot rule out the expanding 
envelope scenario on this information alone because so far we have only considered the 
simplest case of spherical symmetry. Asymmetries can easily arise due to asymmetric 
ejections or structures in the ambient medium. The spectral profiles of the measured clouds 
indicate temperatures of ~ 10 4 K with small turbulence motion. Our freely expanding 
envelope scenario fits these features, too. 

Explanation of low metalicities in these clouds with measured D/H are not a problem 
either, because even though the supermassive stars may undergo hot hydrogen burning 
and break-out to the rp process they may not eject significant amounts of nucleosynthesis 
products from deep in their cores. These stars, however, could produce huge amounts of 
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4 He (FWW). 

Selection criteria dictate that systems with high deuterium abundances are most likely 
to be picked out. Therefore, if Pop III supermassive stars existed in the early stage of our 
universe and they ejected material, all or some of the 3 Lyman limit systems with detected 
deuterium abundances may be accounted for by the ejected envelopes of these stars. The 
required population of these supermassive stars need not be great. One thousandth of 
baryons in the universe forming supermassive stars can be enough. There is no conflict in 
having a universe where <^ 0.1% of its baryonic mass is in ~ 10 5 M Q black holes Some of 
these black hole remnants of supermassive stars may even assume a role in galaxy formation 
by accreting mass and becoming the massive black holes observed at the center of many 
galaxies (van der Marel et al. 1997). 

4. Discussion 

Thus far we have found that: 

• Anti-electron neutrinos from collapse of supermassive stars can produce enough free 
neutrons in outer layers of supermassive stars to boost the deuterium abundance 
there to ~ 1(T 5 to 1(T 4 . 

• If the deuterium-enriched outer layers get ejected, they can reach a large enough size 
that they could mimic low metalicity Lyman limit absorbers at high redshifts. The 
chance of seeing such a Lyman limit absorber in a quasar line of sight is about 1 in 
several tens if 0.01% to 0.1% of baryons are in supermassive stars at a redshift of 2 to 
3, comparable to the current detection probability. 
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Here we especially stress that the ejecta cannot enhance the global deuterium 
abundance in any significant manner, because (1) only a small fraction of baryons may be 
in the form of supermassive stars (otherwise most of the baryons now would be in black 
holes), and (2) only a small fraction of supermassive star material may be ejected due to the 
depth of the gravitational potential. However, we also note that since the chance of seeing 
such a local enrichment of deuterium in quasar lines of sight may not be low, we have to 
be extremely cautious before we can take a D/H measurement in Lyman limit absorbers as 
primordial. So how can we assure ourselves that measurements of D/H in Lyman a systems 
reflect the primordial value? 

First, we do not expect most Lyman limit systems to be envelopes of supermassive 
stars. On the other hand, the enrichment process is sensitive to the conditions in the 
parent stars, and should not yield a uniform deuterium abundance. Consequently, if further 
observations yield significantly more detections with a concordant deuterium abundance, it 
is highly unlikely that we are measuring enriched values. 

Second, if the size of Lyman a absorbers with measured D/H is established by double 
quasar lines of sight to be 3> 1 kpc, they cannot be ejecta of supermassive stars. Their 
deuterium abundances are therefore not enriched. 

Third, the epoch of supermassive stars may not last long. Thus, observations by the 
Hubble Space Telescope of low metalicity absorbers at z < 1 might avoid the enrichment 
problem. However, one should then also worry about the destruction problem, that is, 
whether these absorbers have been processed through stars and their deuterium has been 
destroyed (Jedamzik & Fuller 1997). 

Ultimately, evidence for supermassive star processing through enhanced D/H values 
would suggest that there is some way in which neutrino-exposed ejecta can escape from 
fairly deep in the gravitational potential well of these stars, and the D/H values derived 
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from Lyman limit absorption systems could probe yet another aspect of structure formation 
in the early universe. Along with an enhanced helium abundance, an enhanced deuterium 
abundance in Lyman a absorption systems may signal the existence of supermassive stars 
at high redshifts. 

We thank Scott Buries, Christian Cardall, Karsten Jedamzik, David Schramm and 
David Tytler for helpful discussions. The work is supported by NASA grant NAG5-3062 
and NSF grant PHY95-03384 at UCSD. 
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Figure Caption 



Figure 1. The plot shows the timescale of n(p/y)D and D(p,7) 3 He, and the dynamic 
timescale against the temperature (in units of 10 9 K) in the star, assuming the density 
in the star p oc T 9 3 . The heavier (lighter) lines are for M 5 pa 1(10). The stippled region 
corresponds to thermodynamic and expansion conditions where deuterium can be produced 
and can survive, for M 5 pa 1. The region does not shift much for M 5 pa 10. 



